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We study macroscopic electrical conductivity of a composite made of straight or coiled nanowires 
suspended in poorly conducting medium. We assume that volume fraction of the wires is so large 
that spaces occupied by them overlap, but there is still enough room to distribute wires isotropically. 
We found a wealth of scaling regimes at different ratios of conductivities of the wire, cri and of the 
medium, (J2, length of wires, their persistent length and volume fraction. There are large ranges of 
parameters where macroscopic conductivity is proportional to ((Ticr2)^''^. These results are directly 
applicable to the calculation of the macroscopic diffusion constant of the nonspecific DNA-binding 
proteins in semi-dilute DNA solution. 



I. INTRODUCTION 

Equilibrium and transport properties of composites are 
of great interest because of their importance both in na- 
ture and technology. One usually wants to character- 
ize the composite macroscopically, determining its effec- 
tive properties such as conductivity, dielectric constant, 
magnetic permeability etc., in terms of properties of the 
respective constituents. Most of the theoretical litera- 
ture on this subject dealt with spherical (or single scale) 
inclusionsiiSiiiiSiSii. At the same time, Monte-Carlo sim- 
ulations and experiments revea]i2i2iiS4iiii24iii4 that elon- 
gated, needle or stick like inclusions, can be very effective 
in modifying the properties of materials even at small vol- 
ume fractions. For example, composites made of metallic 
wires with aspect ratio a/l ^ 1 (a is the radius of the 
wire and I is its length) immersed in a good insulator 
exhibit record values of the dielectric constan1>i^. The 
transport properties of such composites were studied in 
many works^iiiiiSii^. A comprehensive review of these 
works and a thorough study of the dielectric response of 
conducting stick composites can be found in Refti^, but 
only in the asymptotic regime of very large conductivity 
of wires. 

Recently, another system attracted a lot of attention. 
It consists of carbon nanotubes dispersed in ceramic or 
plastic material. It was shown that nanotubes can greatly 
enhance the electrical and thermal conductivities of the 
materiaiiiiSiiLima. 

In this paper we present scaling theory of the macro- 
scopic conductivity a of the suspension of well conducting 
wires with conductivity ui in poorly conducting medium 
with a finite conductivity (72 <C cri. The wires can be 
rigid sticks or flexible and coiled like conducting poly- 
mers. We imagine that they are dispersed, randomly 
oriented and frozen in the medium. 

On the first glance, the problem of macroscopic con- 
ductivity of the composite seems to belong to the realm of 
percolation^. Indeed, this would be true for the conduct- 
ing wires in perfectly insulating medium, a2 — 0, where 
macroscopic conductivity could only be realized through 
the direct contacts between wires. In this article, we are 
interested in a different problem - we assume that the 



medium does conduct, (T2 ^ 0, albeit poorly (a2 <C cri). 
In this case, although overall macroscopic current is car- 
ried mostly along the wires, it is still able to switch from 
wire to wire through the medium, depending on the ran- 
dom disordered configuration of the wires. Accordingly, 
we mostly consider volume fraction of wires 4> to be not 
only small (j> <^ 1, but actually smaller than the corre- 
sponding percolation threshold cj) ^ a/ 1, such that the 
direct contacts between wires are rare and completely 
negligible. We show later how our results properly cross 
over to that of percolation at larger (/>. 

For very dilute system of wires, when the distance be- 
tween wires is much larger than the length of the wire, the 
effective conductivity a is trivially close to the conduc- 
tivity of the medium (T2- We therefore mostly deal with 
larger concentrations with > a^/P where the spheres 
containing each wire strongly overlap (see Fig. In 
the parlance of polymer science, we study the semidilute 
systemSi of wires. In terms of increasing concentration, 
our theory continues as long as there remains enough 
room to distribute wires isotropically. We show that al- 
though (/) is small in semi-dilute system of wires, macro- 
scopic conductivity a is dramatically enhanced when 
compared to a2- 

The useful image to think about is a single typical cur- 
rent line in the system. It follows inside one wire for a 
long distance and then bridges to a neighboring wire over 
a more-or-less narrow gap in the medium, and then con- 
tinues again in the wire for a long distance. We denote 
A the length scale of one continuous stretch of current 
line in one wire, it can be called correlation length. This 
is the key concept of the paper. With the increase of 
wire conductivity cti, the correlation length A increases, 
and so does macroscopic conductivity cr, until A gets as 
large as the wire length I, then macroscopic conductivity 
a saturates at values independent on cti. Remarkably, 
for flexible wires there are wide scaling regimes, where 
a oc {(Ji<72Y^^. Such dependence is known for a nar- 
row vicinity of percolation threshold in two dimensional 
isotropic mixturesii^ but, to the best of our knowledge, 
has never been claimed for a broad range of parameters. 

Our theory of effective conductivity of composites can 
be easily applied to a completely different problem, for 
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which the meaning of the correlation length A is particu- 
larly obvious. Namely, we speak of diffusion of proteins 
through the semi-dilute system of dsDNA molecules. 
Many proteins have positively charged domains on their 
surfaces, which provide for nonspecific attraction to the 
negatively charged surface of the double helical DNA. 
Such proteins stick to DNA and diffuse along DNA for 
some time, then get desorbed and wander in 3D, then 
get adsorbed for another tour of ID diffusion, and so 
on. These phenomena are believed to be behind the abil- 
ity of proteins to locate their specific functional targets 
on DNA faster than simple 3D diffusionSSiSiS^. As 
regards macroscopic diffusion of proteins through semi- 
dilute DNA solution, we show later in this article that 
this problem can be easily reduced to that of conductiv- 
ity in the composite with nanowires. As a result, Fig. |21 
part of Figs. 13 and ^ (which is redrawn in the Fig. ]^ 
can all be understood in terms of macroscopic diffusion 
if one uses translation keys provided in the captions of 
these figures. The results are also summarized in Table 
m Clearly, in the case of protein diffusion A is the length 
of one tour of protein diffusion along DNA (including 
episodes of activated desorbtion if they are followed by 
correlated re-adsorbtion). 

Our results are presented by the "phase diagrams" in 
the log-log plane of parameters 4> vs. s = a\l 02 shown in 
Figs. [21iniand2| They specify scaling regimes of different 
power laws formulae for a listed in the Table Rela- 
tively simple phase diagram of Fig. |21 presents results for 
straight wires, while more complicated phase diagrams 
of Figs. 131 and 0] are constructed for semi- flexible wires 
characterized by a large persistence length p, such that 

The plan of this article is as follows. In section HTl we 
first consider the relatively simple case when each wire 
is straight. In this situation we explain the main idea of 
our theory and identify several scaling regimes. We then 
consider a more complicated case when wires are flexible 
and coiled (section IIII|) . We continue in section IIVI by 
using these results for the macroscopic diffusion constant 
of proteins in semidilute DNA system. Finally, we con- 
clude with the discussion of other possible applications 
of our work (section lv)l . 

In this article restrict ourselves to scaling approxima- 
tion for the conductivity and to delineating the corre- 
sponding scaling regimes. In our scaling theory, we drop 
away both all numerical factors and, moreover, also all 
logarithmic factors, which do exist in the problem, be- 
cause it deals with strongly elongated cylinders. 



II. STRAIGHT WIRES 

In this section, we concentrate on the system of 
straight well conducting sticks, such as, e.g., carbon nan- 
otubes, suspended in the medium of lower conductivity 
- see Figure QJi. First of all let us note that percolation 
through the wires and the direct contacts where wires 




(a) (b) 

FIG. 1: Local view of the semidilute solutions of wires. The 
correlation length along the wire is shown in lighter color than 
the rest of the network. The conducting channel is shadowed. 
In a), wires are straight sticks. In b), wires are flexible Gaus- 
sian coils. If the mesh size is not longer than persistence 
length, so the wire within each mesh is essentially straight. 
At lesser density, the mesh size is longer than the persistence 
length, and then the wire in the mesh would be wiggly. 



touch each other starts in such system when volume frac- 
tion of wires exceeds a critical value which is of the order 
oi (p ^ a/l (see Ref»ii*i^). It is not a coincidence that at 
about the same concentration it becomes impossible to 
place wires randomly and isotropicallyS^ - both percola- 
tion and nematic ordering occur at the concentration at 
which there are on average about two contacts per each 
stick. 

As stated in the introduction, we consider the range 
of concentrations /P < (f> < a/l. Only in the very end 
we comment on the role of percolation in our system. 
As long as (j) a/l wires are still oriented isotropically 
and the contacts of wires can still be neglected, i.e., we 
are below percolation threshold; the latter means that 
overall macroscopic conductivity is entirely due to the 
fact that the medium does conduct, (72 7^ 0, albeit maybe 
not too well. On the other hand, the spheres containing 
wires already overlap strongly (0 > a? /P). The latter 
condition means that we deal with semidilute solution 
of sticks - the system that locally looks like a network 
with certain mesh size r (see Figure n^-)- In the scaling 
sense, r is the same as the characteristic radius of density- 
density correlation function, and can be estimated by 
noticing that one stick within one mesh makes density 
about ~ rc? jr^ ~ (/), therefore r ~ a0~^/^. 

Let us start with the simplest case when a\ is not sig- 
nificantly larger than a^- Then current basically has no 
incentive to concentrate into the wires, instead it flows all 
over the place, and effective conductivity is simply that 
of the medium: 

cr ~ (72 • (regime A) (1) 

In the diagram Figure |21 this regime is denoted as A. 

It is also quite easy to find perturbative corrections to 
the result Eq. J^l assuming small "conductivity contrast" 
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{oxIgi) — 1 <C 1. In this case cr ~ (T2 + (fi — o'2)0, which 
follows from the fact that to the first order in perturba- 
tion the current lines remain unaffected by the difference 
between a\ and cr2, they remain parallel straight lines 
(assuming for simplicity the simplest geometry of a uni- 
form current field), and each of them runs through a\ 
material instead of 02 over the fraction of its length. 
This perturbative result suggests that regime A continues 
as long as CTii/) ^ o'2- 

Let us now switch to the more challenging case, when 
CTi is so much larger than 02 that the current is mostly 
carried by the wires. To imagine the flow of current 
through the composite material in this case, it is useful 
to think of a single current line. As stated in the in- 
troduction, such line typically consists of long stretches 
along one wire followed by relatively short switches from 
wire to wire, and the important parameter is the typical 
length over which current line follows inside one wire, we 
denote it as A. 

Our plan is to consider A as a variational parameter. 
That means, we first imagine the current distribution at 
some given value of A, determine the resistance (or con- 
ductance) of the macroscopic sample as a function of A, 
and then try to optimize A accordingly. The justification 
of this procedure comes from the fact that all parts of 
our material obey linear Ohm's law (linear response the- 
ory), so that the requirement of minimal dissipation, or 
maximal overall conductance, is mathematically exactly 
equivalent to the Kirchoff 's laws determining distribution 
of currents in the network of resistorg^^. 

To begin with, let us consider the most interesting case 
when A is much shorter than total length of one wire, 
but still larger than the mesh size r: r <C A ^ Z. Con- 
sider then a cube of the size about A inside our macro- 
scopic sample. On the one hand, overall conductivity 
on the scale of this cube is already about the same as 
that of a macroscopic body, we denote it cr; the resis- 
tance of one A-size cube is, in other words, about 1/Acr. 
On the other hand, we can estimate this resistance con- 
sidering wires inside the cube. There are about }? jXr^ 
of the wires crossing the cube, because the distance be- 
tween wires is about r and, therefore, each wire can be 
thought of as dressed in a sleeve of the thickness about 
r and volume about Ar^ (see Fig. ^). The sleeve can 
be thought of as a weakly leaking insulation for the wire. 
Each piece of the sleeve of the length about r bridges 
given wire to another one through the resistance about 
r/u2r'^, and about A/r such bridges are connected in 
parallel, yielding overall resistance connecting the given 
wire as (l/(T2r)(r/A) — 1/(T2A. This is connected in se- 
ries with the wire itself, producing a conducting channel 
with resistance X/a^ai + \/cr2X- Since all (or sizeable 
fraction of all) A^/r^ conducting channels in the cube 
are in parallel, we finally arrive at the cube resistance 
as (r^/A^) [X/a^ai -f 1/(T2A). Equating this to l/Acr, we 
arrive at the following estimate of effective macroscopic 



conductivity 



II ai + a2/^2A2 • ^ > 

Let us now analyze this result. As a function of A, 
macroscopic conductivity does not appear to have max- 
imum at any finite A. This does not necessarily mean 
that A is going to increase ad infinum: rather, it sug- 
gests that more accurate calculation is required to de- 
termine A. Luckily, such more accurate calculation is 
not necessary to determine the quantity of our interest - 
macroscopic conductivity a. Indeed, as soon as A exceeds 
certain threshold, namely 

A > ail 02 , (3) 

the second term in denominator of formula , which de- 
scribes resistance of wire-to-wire bridges, becomes sub- 
dominant, and must be neglected within the accuracy of 
our scaling estimates. This yields 

(T ~ cri(/) . (regime B) (4) 

As expected, regimes A and B cross over smoothly on 
the line cti/(T2 ~ 0- 

In the regime B, macroscopic conductivity does not 
depend on CT2, conductivity of the medium. This hap- 
pens because (J\Ig2 is so large that current mostly fiows 
through the wires, but at the same time g\ is not large 
enough to make resistance of wires insignificant and thus 
unmask the resistances of narrow gaps between wires. 
This is also why the resulting macroscopic conductiv- 
ity is in this regime insensitive to the exact value of A: 
overall distance travelled by any particular line of cur- 
rent through the wires scales simply as the sample size 
and, to the scaling accuracy, does not depend on how 
frequently current switches from wire to wire - precisely 
because wires are straight. 

To complete our analysis of the regime B, let us note 
that the quantity o\l 02 which appears in the formula 
can be understood in the following way. Imagine 
one straight wire in an infinite medium, and suppose we 
somehow feed current into a certain point of this wire. 
Current will flow away from 0, mostly along the wire, 
but also slightly leaking into the environment. Due to 
this leaking the current remaining in the wire will decay 
exponentially as we move away from 0, with the decay 
length equal to ay/oijo^. Indeed, the length of current 
decay for a single wire is estimated by the condition that 
the resistance of wire over the length A, which is about 
XIgio?, should be about the same as the resistance of 
the medium "in perpendicular direction" , which is about 
I/CT2A; equating these two returns the result P|). 

Let us now switch to the next regime H, which arises 
because of the finite length of sticks. Indeed, the length A, 
over which current flows in one wire, cannot exceed the 
wire length I. According to formula this crossover 
happens along the border s = PI a?. Indeed, starting 
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from this s, we get even a a/ u\I > I, and so A cannot 
satisfy Eq. Q. To find conductivity in the regime H, we 
should replace A by Z in Eq. Since the second term 
in the denominator dominates, we arrive at 

a ^ (I){l/afa2 ■ (regime H) (5) 

We get the same P dependence as predicted in 
referencei^. This effective conductivity has no depen- 
dence on ai because cti is so high that overall resistance 
is entirely concentrated in the narrow gaps where current 
should switch from wire to wire. As a result, cti does not 
enter the formula. 

If we decrease (j) and look at the dilute regime (j) < 
a? /P, our scaling theory suggests that cr ~ cr2- More 
accurate analysis of this regime was performed in the 
work^^ using first order perturbation theory in powers 
of (/), which is applicable in the dilute regime, as long as 
current field around one wire does not affect the neigh- 
boring wires. The achievement of that work is that the 
authors were able to take into account, to the first order 
in their perturbation theory, both numerical coefficients 
and logarithmic factors proportional to ln(//a) (in our 
notations). Up to these factors, which we systematically 
neglect, all our scaling results cross-over smoothly with 
the results of the work^*. Moreover, the simplified ex- 
pressions of perturbation results are given in the work^® 
in the form of three formulae, describing the dilute so- 
lution in the ranges (in our notations) / ^ a^ai/<j2, 

I K, a^J f7i/(T2, and I <C a\J (Ji/cr2\ this sheds an ad- 
ditional light on the importance of the lengths scale 
a\J a\l 02 introduced by us here in the present work. 

Thus, we have completed our consideration of the 
phase diagram up to the concentration about </) = ajl. It 
would be frustrating to stop at this point, because, for ex- 
ample, the experiments with suspensions of carbon nan- 
otubes often use loadings with > ajl to achieve larger 
electrical and thermal conductivities. Let us therefore 
discuss what happens to conductivity if one can man- 
age to create isotropic suspension with </> > ajl. Sup- 
pose the length of wire between direct contacts, or the 
mesh size of the percolating network, is C- Then, the 
number of electrically parallel wires within a cube of size 
C is of the order of C'^(/)/(Ca^), each with the resistance 
about C/(cria^). Therefore the total resistance of the 
cube scales as l/(criC'^), yielding the effective conductiv- 
ity about (710. In the scaling sense, this is the same a 
as that in Regime B. Thus, regime B continues to higher 
concentrations (\) > a/ 1. Trivially, the lower bound of this 
continued regime B remains to be the condition s > 1/0, 
since at lower s current does not concentrate in wires, 
percolation gives no help. More interestingly, percolation 
produces no effect on conductivity a as long as s < (l/a)^ 
because the percolating network gives no advantage over 
the conducting channels made of combination of wires 
and surrounding medium. Percolation does have effect 
when s > (Z/a)^, where conductivity as a function of 
concentration (j) changes rapidly from a ~ 0(Z^/a^)cr2 
below percolation to cr '--^ (fxri above. This change occurs 
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FIG. 2: Diagram of scaling regimes for the case of straight 
wires (e.g. carbon nanotubes). Both volume fraction of the 
wires (j) and s = axja-i a:xes are in the logarithmic scale. Each 
line on this plane marks a cross-over between scaling regimes. 
The widened line at </!> = ajl shows the location of conduc- 
tivity jump around the percolation threshold. Summary of 
the conductivities is provided in Table If s is replaced by 
yd, the diagram represents macroscopic diffusion constant of 
protein discussed in section HVI 



around the threshold ip = (pc — a/l^ over the interval 
of width of order (pc, which is schematically plotted as 
the shaded cone in the diagram Fig. [21 In this range, 
the conductivity has critical behavior similar to that dis- 
cussed in Ref.^. The detailed structure of this range is 
beyond the scope of this paper. 

We should emphasize that in all our considerations 
we completely disregard the resistivity of the contacts, 
either between a wire and surrounding medium, or be- 
tween two wires in contact. In particular, everything we 
said about percolation assumes that whenever there is a 
touch between two wires, it presents an electrical con- 
tact of vanishing resistance. In fact, this assumption is 
model sensitive and it is not necessarily good in practi- 
cal cases. For instance, our theory predicts that above 
percolation threshold, when p > a/l, effective electric 
conductivity cr grows linearly with (p, independently of 
s being larger or smaller than {l/a)'^. Similar predic- 
tion holds also for thermal conductivity. However it is 
not compatible with the apparently super-linear growth 
of effective electrical or thermal conductivity observed in 
experimentsiiiSiiSiSS. It is not clear whether interfacial 
resistance can help to explain these experimental data. 
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TABLE I: The summary of macroscopic conductivities and 
diffusion constants in various regimes. 
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FIG. 3: Diagram of scaling regimes for the case of flexible 
Gaussian coiled wires (e.g. conducting polymers) with length 
p < I < p'^/a. Summary of the conductivities is provided in 
Table n 




FIG. 4: Diagram of scaling regimes for the case of flexible 
Gaussian coiled wires with length p'^/a < I < p? ja^ . 



III. GAUSSIAN COILED WIRES 

When wires are flexible (e.g. conducting polymer), our 
theory developed in section |n] needs modifications to ac- 
count for different fractal properties of the wires. We con- 
sider the system of semi-flexible conducting polymers as 
an example (see Fig. QJ)) . We assume that each polymer 
has contour length /, relatively large persistence length 
p and the radius (thickness) a. We assume I ^ p ^ a. 
Throughout this work we require I < /o? and this lets 
us disregard the effect of excluded volume on the polymer 
statistics, considering polymer coil as Gaussian. There- 
fore, if we take a piece of polymer of the contour length 
A, then its size in space scales as 



A when \ < p 
^/\p when \ > p 



(6) 



The overlap of coils starts when the volume fraction of 
polymers exceeds volume fraction of the chain inside one 
coil: - la'^/{lpf/'^ = a? / {I'^/'^p^/'^). The nematic or- 
dering of wires starts at a larger volume fraction </> ~ a /p. 
The percolation threshold can be estimated from follow- 
ing argument. Percolation happens when each wire has 
roughly two direct contacts with other wires. Following 
RefiSi, the number n of contacts per wire is of the order 
of 



n ^ (fil/a 



(7) 



(each piece of the polymer of the length a has probability 
(j) to be in touch). Requiring that n is of the order of 1, we 
obtain the percolation threshold (j)c ~ a/ 1. Percolation 
starts earlier than the nematic ordering and 4>c, because 
percolation requires a couple of contacts per wire, while 
nematic ordering requires a contact per a smaller length 
of a statistical segment {'^ p). Therefore, 0c divides into 
two parts the range of concentrations when the system 
is semi-dilute but isotropic j ip-^'^p^l'^) < < ajp: 
^2^^-^1/2^3/2-1 < < 0^ where we neglect the effect of 
direct contacts; and (f)c < <(> < a/p where the effect of the 
percolating wires must be included in our theory. 

Let us start with determining the mesh size 7i2i. Sup- 
pose that the polymer within each mesh has a contour 
length g. ft makes a density about ~ go? /r^ which 
must be about the overall average density 0. Thus, 
go} jr^ ~ 0. Second relation between g and r depends on 
whether mesh size is bigger or smaller than persistence 
length p: 



g \i g<p 
V9P ''^i 9 > P 



(8) 



Accordingly, one obtains 



9 ' 
9 ' 



, a^^ if 
°' , r ^ — if 



- > d>> ^ 



(Ip)^/' 



< < 



(9) 
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The upper hue corresponds to such a dense network 
that every mesh is shorter than persistence length and 
polymer is essentially straight within each mesh (see Fig. 
^). The lower line describes much less concentrated 
network, in which every mesh is represented by a little 
Gaussian coil. Depending on the relation between I and 
p, percolation can start before or after (j) ~ a^/p^. Our 
results for these two cases are summarized in Figs. |21 01 
and in the Table |l| 

When the correlation length X < p, the polymer within 
correlation length is straight. So we can directly apply 
what we got for the straight wire case and obtain regimes 
A and B. 

For other regimes with X > p, the derivation has to be 
performed from the beginning. So, we consider a typi- 
cal cube inside the composite with size ^ such that every 
polymer enclosed in this cube has contour length of the 
order of correlation length A. There are about r^^xfg-j 
electrically parallel conducting channels in this cube (be- 
cause X/g is the number of meshes visited by one wire, 
and is the volume of each such mesh). Each chan- 
nel consists of the wire itself, with resistance A/(cria^), 
and the wire is connected in series with the group of X/g 
parallel connected bridges, each of resistance r/{<T2r'^)- 
Thus, the resistance of the cube scales as 



X/{aia^)+g/{a2Xr) 
5CV(Ar3) 



(10) 



which should be equated to 1/ (a^) . Therefore we obtain 

a^g/r^ 



Xy{a,e)+9ayia2re) 



(11) 



Once again, if the wire remains straight over the length 
A, so that ^ ~ A, we are back at the situation described 
by formula |2Jl, and we can reproduce the corresponding 
result for the regime B, Q. More interestingly, we can 
now consider the case when A > p and the polymer of 
the length A is the Gaussian coil: ^ ^ In this case 

we obtain 



pa^g/r^ 



A/<Ti + a'^g/ ((T2rA) 



(12) 



Now conductivity has the well defined maximum at the 
well defined value of A: 



(72 r 



1/2 



(13) 



Not coincidentally, this result for the correlation length 
A corresponds to equating two terms in the denomina- 
tor of Eq. 1)12(1 - the resistance of wire with correlation 
length A and the resistance of the surrounding "sleeve" 
of thickness r in the medium. 

Plugging A from Eq. ifT^ back into Eq. (fT^ and 
applying Eqs. (Q, we obtain the following two scaling 
regimes: 



Regime C, where polymer on the scale A is Gaussian 
(lower line in the Eq. lO, but the polymer within each 
mesh is still straight; 

Regime D, where polymer is Gaussian even within each 
mesh (lower lines in the Eqs. ©). 

When A reaches the entire length of the polymer I, 
we should replace A by Z in Eq. H12|) . Then the second 
term in the denominator dominates. Since we have two 
different kinds of meshes represented by formulae , we 
have two more regimes: 

Regime E, where A = I and the polymer within each 
mesh is straight (upper lines in the Eq. Ipjjl): 

Regime F, where A = Z, but the polymer within each 
mesh is Gaussian (lower lines in the Eq. (|5J). 

We should emphasize that regime E exists only in the 
case where / < p^/a (see Fig. OJ. If I > p'^ /a, percolation 
starts so early that the polymer within each mesh is a 
Gaussian coil. This case is plotted in Fig. 0] 

When we increase ai/a2, the effective conductivity 
grows from tT2 (regime A) and finally saturates in regimes 
E and F with values having no dependence on ai. As 
we discussed, for regime A, transport through wires is 
not at play while in regimes E and F, di is so large 
compared to (72 that the wires are effectively "super- 
conducting". More interestingly, for broad ranges of (jj 
and (7i /(72 (regimes C and D) , cr is proportional to ^Ja\a2. 
Such dependence is known for a narrow vicinity of per- 
colation threshold in isotropic mixtures^ but, as far as 
we know, it has never been noticed for a broad range of 
parameters. The width of the range grows as Z^. 

When the volume fraction is larger than the percola- 
tion threshold, the effect of percolating wires can not be 
neglected. If (Ji/a^ is relatively small, the transport of 
current is mainly realized through the conducting chan- 
nel we have discussed. But when oxja^ is large enough, 
percolation through the directly connected wires domi- 
nates. The crossover is determined by equating the con- 
ductivity due to the untouched wires and surrounding 
medium and the conductivity due to percolating wires. 
We have already calculated the first conductivity. The 
later one can be calculated by the following argument. 
Let us denote the length of wire between contacts by 
C,. It can be estimated as C '^-^ l/{4'l4'c) ^ a/0- Be- 
cause we require (j) < a/p, C, is larger than the persistence 
length p and thus the distance it covers in space scales 
as ~ (Cp)^^^- Within a cube with size (Cp)^^^j there are 
(^p)^/^(/)/(^a^) ~ (p/a)^/^(/)^/^ electrically parallel wires. 
So the conductance scales as aipa/Q. It can be also ex- 
pressed using effective conductivity, it is a {CpY^^. Com- 
paring these two conductances, we obtain a {p/a)(j>'^ai, 
which is the effective conductivity in regime G. It crosses 
over smoothly to the regimes E and F (E only exists for 
case I < p'^, which is represented in Fig. O. One can 
also obtain the border by equating the correlation length 
A to the length between direct contacts C,. Since we as- 
sume the current can switch wires freely at the contacts, 
A can not grow above C. On the other hand in the scal- 
ing approach we use, there is a conductivity jump around 
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FIG. 5: Diagram of scaling regimes for the case of protein 
diffusion in Gaussian coiled DNA solution. Summary of the 
macroscopic diffusion constants is provided in Table H] 



the percolation threshold 0c ~ o/^, which is plotted as 
the widened line in both figures. Actually the jump of 
conductivity is eliminated when we consider the critical 
behavior of conductivity at — 0c ^ 0c but our scaling 
theory is not designed to see such details. 



IV. MACROSCOPIC DIFFUSION CONSTANT 
OF PROTEINS IN SEMIDILUTE DNA SYSTEM 

The theory we have developed for the effective con- 
ductivity of composites can be used to study the macro- 
scopic diffusion constant of the nonspecific DNA-binding 
proteins in semidilute DNA solutions. 

For simplicity, we make the following assumptions: (i) 
protein can be non-specifically adsorbed on any place of 
DNA; (ii) non-specific adsorption energy e, or the corre- 
sponding constant y — e'^l^^'^ ^ are the same everywhere 
along DNA (sequence-independent); (iii) every protein 
molecule has just one surface patch capable to stick to 
the DNA, so proteins do not serve as cross-linkers for 
the DNA; (iv) non-specifically bound protein can diffuse 
along DNA with the diffusion coefficient while pro- 
tein dissolved in surrounding water diffuses in 3D with 
diffusion constant (v) while protein is diffusing, the 
DNA remains immobile. 

To make the dictionary of translation between con- 
ductivity and diffusion languages, the easiest way is to 
step up the generality in writing down the expressions 
for current density j in either conductivity or diffusion 
problem. In both cases, as long as we consider stationary 
process, current is subject to the no-divergence condition: 
divj = 0. For the electric current driven by the potential 



gradient. Ohm's law reads j = — (t(x)V0; and for the the 
diffusion problem, current driven by the gradient of total 
chemical potential is described by similar Smoluchowski 
equation: j — — D(x)c(x)V (Inc(x) — e(x)/fcBT). Here, 
we assume for a moment, that electrical conductivity, dif- 
fusion coefficient, protein concentration c, and binding 
energy e or y = e!^I^BT j^g^yg ^ some general dependence 
of space coordinates x. In fact, for our case, this space 
dependence is very simple: within narrow regions along 
the wires or along DNA, we have ct(x) = a\ , and similarly 
D(x) = D\ and ?/(x) = y (remember that y — e^/*^^^); 
for all other places x we have <t(x) = (T2, £'(x) = D3 
and j/(x) = 1. As regards concentration, diffusion equa- 
tion also implies (since chemical potential is continuous) 
that locally there is the equilibrium relation between the 
ID concentration of non-specifically adsorbed proteins, 
ci, and concentration of proteins remaining free in the 
nearby solution C3: 



c\l[cza^) = y 



(14) 



where a is the length scale such that ci/a^ is the 3D 
concentration of proteins within the region around DNA 
where proteins are adsorbed. Comparing the equations, 
we see that complete mapping is achieved by the substi- 
tutions (j\ <-> DiCi/ a?' , a-i <-!■ -D3C3. Similarly writing the 
effective macroscopic equations in terms of macroscopic 
conductivity a and macroscopic diffusion coefficient Z?, 
one finds a <-^- Dcz- In terms of more convenient di- 
mensionless quantities, and taking into account the local 
adsorbtion equilibrium (|14|l . these rules read: 



(Ti Dx a D 

(72 D3 ' (72 -D3 



(15) 



We can, therefore, directly address macroscopic diffusion 
based on our results for macroscopic conductivity. Sub- 
stituting Eq. H15() into our results for a, we obtain the 
macroscopic diffusion constants of the proteins in the 
DNA solution expressed in terms of Di and D3 for all 
the regimes except regime G, for which the applicability 
of percolation results to the DNA case is doubtful.^" 

Thus, Fig. 121 can be used for the macroscopic diffusion 
constants for the straight DNA case if we replace ai/(72 
by yDi/D^. However, for the gaussian coiled DNA case, 
without regime G, Figs. 13 and 0] should be modified. 
Resulting phase diagram is shown in Fig. |S1 One can 
easily get this figure from Fig. Olby removing the borders 
the regime G makes with other regimes and extending 
the border line between regimes E and C to the right 
boundary of the phase diagram. The results for various 
scaling regimes are also summarized in the right most 
column of Table 

Measuring macroscopic diffusion of proteins is a 
promising way to test our predictions, ft is therefore 
useful to comment a little deeper on the nature of macro- 
scopic diffusion coefficient D. The way it is defined above 
is adequate for a macroscopic experiment, because D es- 
tablishes the proportionality between the flow of proteins 
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and the gradient of concentration of dissolved proteins. 
In such experiment, the presence of a large number of 
proteins adsorbed on DNA is not directly relevant. How- 
ever, in a different experiment, for instance, in tracking 
the random walks of a single protein molecule, a differ- 
ent diffusion coefficient, D, will be relevant, such that 
Dc — Dc3, where c is total concentration of proteins, in- 
cluding the adsorbed ones: c = (ci/a^)0-|- 03(1 — (/)). Us- 
ing (dl), one then gets D ^ D/ {1~ + y(f>). The differ- 
ence between these two diffusion coefficients is marginal 
when absorbtion is weak (y <C 1 and <C 1), but it 
becomes very much pronounced when the absorbtion is 
strong, or y is large: D ~ D/{y(j)) <C D. This result has 
simple physical meaning: every particular protein will be 
adsorbed most of the time, so its diffusive motion will be 
slow, but overall flow of proteins will not be that slow 
because of a large number of proteins. 

Using the macroscopic diffusion constant, we can also 
re-derive the results of the work^^ concerning the rates 
of protein searching for specific places on globular DNA. 
Thus, measuring D or D is another way to verify the 
results of Ref.^ in those crowded regimes. 

V. CONCLUSION 

In this paper we studied a plethora of different scal- 
ing regimes for conductivity of a suspension of wires in 
poorly conducting medium. Our results are applicable 
to suspensions of metallic wires in poorly conducting 
medium at room temperature. In this case our generic 
description of the system only by two macroscopic local 
conductivities is valid, because typically the surfaces of 
nanowires are so dirty that any surface barrier for elec- 
trons of the metal is sufficiently well conducting due to 
hopping through localized states. 



We also mention carbon nanotube suspensions as a 
possible application of our theory. In this case one may 
worry about the role of the surface resistance on the 
nanotube-medium interface, so that our theory strictly 
speaking only estimates the effective conductivity from 
above. Including surface resistance or allowing for influ- 
ence of environment on conductivity of nanotubes would 
require new parameters in the theory, making it much 
more complicated. We believe that both real and com- 
puter experiments should be first compared with the sim- 
plest and generic model presented here in this paper be- 
fore one starts developing more complicated theories. 

The serious advantage of our generic theory is that 
it can be applied in a variety of different problems be- 
yond electric conductivity, for instance, to thermal con- 
ductivity of well thermally conducting wires in a weaker 
thermally conducting medium, to macroscopic dielectric 
constant of suspended metallic wires and to wires with 
large magnetic susceptibility. We also applied our theory 
to the calculation of the macroscopic diffusion constant 
of the nonspecific DNA-binding proteins in semi-dilute 
DNA solution. 

The latter application is also promising in terms of 
computational tests of our theory along the lines of the 
recent workS^i^S. 
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As we mentioned, even for the wires the idea of direct con- 
tact is very much model-dependent. It is even more so for 
the protein and DNA case, because in this case "contact 
between wires" should mean the possibility for the protein 
to switch from one DNA to the other without activation. It 
might be possible in some systems but impossible in oth- 
ers; besides, there are quite a few other effects which are 
beyond our theory, such as excluded volume constraints 
for the proteins which becomes significant when two DNA 
pieces are close by - protein may have difficulties difi^using 
along one of them, like a big truck under a low bridge on a 
highway. We do not consider all these questions in this pa- 
per, and only consider DNA system well below percolation 
threshold. 



